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Abstract 

We study the sliding state of a two-dimensional Wigner crystal in a strong 
magnetic field and a random impurity potential. Using a high-velocity per- 



^r , turbation theory, we compute the nonlinear conductivity, various correlation 

functions, and the interference effects arising in combined AC + DC elec- 
tric effects, including the Shapiro anomaly and the linear response to an AC 

field. Disorder is found to induce mainly transverse distortions in the slid- 
es ■ 

ing state of the lattice. The Hall resistivity retains its classical value. We 

S^ . find that, within the large velocity perturbation theory, free carriers which 

affect the longitudinal phonon modes of the Wigner crystal do not change the 

form of the nonlinear conductivity. We compare the present sliding Wigner 

crystal in a strong magnetic field to the conventional sliding charge-density 

wave systems. Our result for the nonlinear conductivity agrees well with the 

I — V characteristics measured in some experiments at low temperatures or 

large depinning fields, for the insulating phases near filling factor v = 1/5. 

We summarize the available experimental data, and point out the differences 



among them. 
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I. INTRODUCTION 

The T = phase diagram of a two-dimensional (2D) electron gas in a strong magnetic 
field is very rich. Two of the possible competing phases are the Wigner crystal phase, which 
is believed to occur at sufficiently strong magnetic field or low density, and the fractionally 
quantized Hall effect (FQHE) liquid phase, which occurs if the field is not too strong and 
the filling factor v has an odd denominator (T|Q|. The two-dimensional electron gas has been 
studied experimentally in modulation doped semiconductor heteroj unctions. Currently, the 
FQHE state is firmly established down to filling factor v = 1/5 || for samples with n-type 
dopants and down to v = 1/3 for samples with p-type dopants. For v slightly greater than 
1/5, but less than 2/9, and for v smaller than 1/5, insulating phases have been observed in 
the cleanest n-type samples currently available 0-0. There is also increasing experimental 
evidence that similar reentrant insulating phases appear in p-type samples for v slightly 
greater and for v less than 1/3 ||. Most of the experimental data on the insulating phases 
in this regime are interpreted as due to the Wigner crystallization 0-0. 

While the insulating nature of the states around v = 1/5 has been clearly demonstrated, 
it remains controversial if the insulating behavior is primarily due to interaction-induced 
Wigner crystallization or primarily due to disorder. Nonlinearities in the conductivity have 
been argued by several groups [|]-[7[] to imply that the insulating phases are due to the 
formation of the Wigner crystal which is then pinned by disorder. However, there exist some 
puzzling discrepancies among experiments under seemingly similar experimental conditions. 
For example, the apparent depinning field (see Eq. 3) differs by almost three orders of 
magnitude between two experiments 0J| on samples with similar zero field mobility. In 
another experiment J7| , nonlinearity in the conductivity was observed in both of these field 
ranges; but the nonlinearity at larger fields was attributed to electron heating effects at 
larger currents. Even within the small field (~ 0.2 mV across samples with typical size ~ 
3 mm) nonlinear conductivity studies, some significant differences exist among experiments 
HJBIJTI]. We shall discuss these and other related experiments in more detail later in the 



paper. 

In view of the experimental situation, it is desirable to provide a theoretical framework 



for describing the properties of a sliding Wigner crystal in a strong magnetic field ]T0| . The 
principle theoretical issue is the proper treatment of the disorder which distorts the Wigner 
crystal during sliding. In this paper, we apply techniques originally developed for studying 
sliding charge-density wave (CDW) systems HTTHl5|] to the present system. We use an elastic 
model for pinning and for the interaction of the sliding Wigner crystal with disorder, similar 
to that used in the Fukuyama- Lee- Rice (FLR) model for CDWs fl6|], and the high- velocity 



perturbation theory [|T2| , p~3|| to calculate various properties of a sliding Wigner crystal in the 
presence of a strong magnetic field and disorder potential. This approach is expected to be 
valid as long as the sliding velocity is large, in comparison to disorder effects, but not so 
large as to invalidate the elastic medium theory. In the context of CDW systems, similar 
perturbation theory [12||13| was compared to numerical simulations, and was found to be 



quantitatively reliable for fields greater than ~ 2 times the threshold field. 

The physics of a Wigner crystal in a strong magnetic field differs that of a conventional 
CDW because: 1) the presence of a large magnetic field gives rise to a large Hall component of 
the motion; 2) the motion is two-dimensional whereas it is one-dimensional in a conventional 
CDW system in arbitrary spatial dimensions, and the related issue of the shear modulus of 
a Wigner crystal being much smaller than the bulk modulus; 3) the presence of a weakly 
screened, two-dimensional long-range Coulomb interaction; 4) the order-parameter, i.e., the 
Fourier-transformed charge density, is non-zero on a two dimensional grid, and does not 
decay very fast with increasing reciprocal lattice vector G. The first point does not affect 
the static properties of the Wigner crystal. Points 2 and 3 have been shown to affect the 
static properties, and together with point 1, the fc-dependence of the collective excitations 
of a Wigner crystal in a significant way [XT] . In the sliding regime where we are interested 
in the dynamic response of the Wigner crystal, they will all affect the results. 

Some recent experiments done at relatively high temperatures J7J have also indicated 
the possible involvement of the thermally excited free carriers in the sliding Wigner crystal. 



For example, the differential conductivity continues to be thermally activated above the 
sliding threshold in Ref. 0, up to a second threshold field, thus raising the possibility that 
the thermally excited free carriers play an important role in determining the conductivity 
of the sliding state. We have therefore considered, within the same perturbation theory 
framework, the effect of coupling to a sliding Wigner crystal a set of free carriers which 
are phenomenologically characterized by a conductivity tensor. We find that the present 
perturbation theory is unable to account for the experimental results in Ref. J7| between the 
two threshold fields. This suggests that the starting point of our theory, i.e., the assumption 
that we can perturb around a uniformly sliding state, is probably not valid in this regime. 

We have also considered the effects of combined AC + DC fields, where interference 
between the internal "washboard" frequency of the sliding motion and the AC driving fre- 
quency can be used to further identify the nature of the sliding state of the insulating phases 
observed experimentally. These include the Shapiro anomaly in the DC conductivity in the 
presence of an AC field, and the AC response in the presence of a large DC current. These 
effects are also studied using the same perturbation theory, carried out to the same order in 
disorder potential as that for the DC nonlinear conductivity. 

It has also been claimed recently that a pinned Wigner crystal in zero external magnetic 
field exists at Si/ ' S1O2 interfaces |H|]. For completeness and convenience of reference, we 
also give our results for zero magnetic field, although they are similar to those for the flux 
lattice sliding motion in a type-II superconductor as studied by Schmid and Hauger and 
by Larkin and Ovchinnicov in Ref. |15[ . The main difference here is that for the present 
Wigner crystal case the long-range Coulomb interaction yields a fc^-dispersive longitudinal 
phonon mode, and this dispersion could be affected by the presence of free-carriers in ways 
that are described in |I7| and also discussed later in this paper. As we will show, however, 



the main distortions in a sliding Wigner crystal are transverse; and many properties of a 
sliding Wigner crystal are not affected by the longitudinal mode spectrum. 

The balance of the paper is as follows. In Sec. II, we introduce our model, make some 
general remarks regarding the Hall resistivity, and summarize recent experiments concern- 



ing the Hall resistivity for the insulating phases around v = 1/5. In Sec. Ill, we study the 
nonlinear conductivity of the Wigner crystal in the sliding regime in the absence of free car- 
riers, treating the scattering by impurities as a perturbation. Various correlation functions 
of the sliding motion are studied in Sec. IV. The effects of the disparity of the bulk modulus 
and the shear modulus are illustrated. In Sec. V, we consider the effects of free carriers 
on the sliding dynamics of the Wigner crystal. In Sec. VI, we consider the two AC + DC 
interference effects mentioned above. Results for the zero magnetic field case are given in 
Sec. VII. All of our main results for the transport properties of a sliding Wigner crystal 
are summarized in Sec. VIII. In Sec. IX, we discuss in some detail the current experimental 
situation and attempt to make contact between our theory and these experiments. Sec. X is 
the conclusion. Readers who are only interested in the results and the current experimental 
status may go to Sees. VIII, IX and X directly. 

II. EQUATION OF MOTION AND THE HALL RESISTIVITY 

In Sec. IIA, we first establish the equation of motion. We remark briefly on its similarities 
to and differences from the equations of motion of related models. In Sec. IIB, we give results 
for the Hall resistivity that can be obtained immediately by examination of the equation 
of motion. The steady state solution to the equation of motion in the absence of random 
potential scattering is discussed in Sec. IIC. 

A. Model and Equation of Motion 

We begin by considering the various forces acting on an element of the Wigner crystal 
at (r, t) in two dimensions, which we treat classically as an elastic medium. We denote the 
two-dimensional displacement of the Wigner crystal at (r,t) as u(f,t), average mass and 
charge densities as p m and p c , lattice constant of the Wigner crystal as a, and background 
dielectric constant e. The following terms enter the equation of motion: 

1) inert ial force: p m u; 
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2) dissipation: Xp m u. A is a phenomenological constant describing damping due to 
everything other than the disorder potential that is treated explicitly below; 

3) Lorentz force: p m u c e z x u, where u c = p c B/p m is the cyclotron frequency, e z is the 
unit vector in the direction of the perpenticular magnetic field; 

4) elastic restoring force: / D(r — r') ■ (u(r', t) — vt)d 2 r'. Here D is the real-space dynamic 
matrix tensor, and vt is the uniform component of u(f, t) in the sliding state of the Wigner 
crystal defined as vt = limT^^limy^^ J ^ / *«(r, t). Upon Fourier transformation, the 
component transverse to k is T) T (k)/p m = c T k 2 and the longitudinal component Y) L (k)/p m = 
c L k. These are valid in the long wavelength limit, i.e., for \u(r)—u(r')\ much less than \r— r'\. 
Following Ref. [ 17fl , we have: c T = Q 2 a 2 a, and the fc-dependent c L = Q 2 (2tcci + (1 + a)a 2 k), 



where a ~ 0.02 is the ratio of the shear modulus with the bulk modulus, and Q 2 = p 2 /p m ae is 
a characteristic frequency of the Wigner crystal, on the order of the zone-boundary phonon 
frequency at zero magnetic field. The first term in c L is due to the long-range Coulomb 
interaction. Here we view the Wigner crystal as a deformable classical solid characterized 
by its dynamic matrix D(A;), given by, for example, Bonsall and Maradudin |Tj|. Quantum 



mechanical corrections have been studied within the time-dependent Hartree-Fock theory 
and are found to be on the order of 10% for the long wavelength phonons that we are 
interested in [ p0|1 ; 

5) disorder potential which couples to the r-dependent charge density p(r): 
p{r)V r Vdis{r + u(f,t)). Here Vdi S (r) is the disorder potential at f. Only the correlation 
function of Vdis will enter our results, so we will not specify its form. In cases where Vdi S is 
due to the remote dopants at a distance d away from the plane of the 2D electron gas, the 
Fourier transformation of V^ s (r), Vdi S (q) decays as e~ qd from simple eletrostatic considera- 
tions; 

6) and finally, the external driving force on the total charge density p c E ext . 
Combining the six terms, we obtain the equation of motion in real space: 



p m u + \p m u + p m uj c e z xu+ D(f - ?) ■ u(r',t)d 2 r' = p c E ext - p(r)V 'fV dis {r '+ u(f,t)). 

(1) 

In a perfectly ordered Wigner crystal, the charge density is determined by its Fourier 
components p{G) at the reciprocal lattice vectors G, via, p{f) = J2q p(G)e lG ' r . In the 



systems of interest, the magnetic field is so strong that the magnetic length Ib = Jhc/eH 
is smaller than the lattice constant, and p(G) decays as p(G) = e~ G 1 bI 2 , so many reciprocal 
lattice vectors will be important, in contrast to the CDW case in which only the lowest G 
is important. 

For later reference, we write down here also the equation of motion in the FLR model 
for a CDW. Apart from the absence of a magnetic field, the most important difference is 
that because of the quasi-one dimensional crystal structure of CDW systems, the electrons 
are free to move only in one direction, which we take to be z. Therefore: 

/-> -» -> 

D zz (f- r')u z (r',t)dr' = p c E e z xt - p(r)—V dis (r + e z u z (r,t)). (2) 

Dimensionality, which is crucial in determining the functional form of the nonlinear conduc- 
tivity of a sliding CDW |H|, enters through the elastic term in the above equation. 



In the absence of the external potential, very weak disorder can pin the overall phase of 



the charge-density wave [|16| . In this case, the phase fluctuations diverge beyond a length 



£ T which defines a typical domain size in a pinned CDW. Similar considerations apply 



here to the Wigner crystal case as well [17]. As a result of pinning, the linear response 
conductivity vanishes. But when E ext is greater than some threshold field Exh-, the pinned 
Wigner crystal will begin to slide. The pinning-depinning transition between the static state 
and the sliding state in a CDW system has been argued to be second order, exhibiting the 



critical phenomena 21 



In the weak pinning limit, following the Fukuyama- Lee- Rice argument |TB,T7], one can 
relate E^h to the static correlation length £ T of the pinned Wigner crystal as follows. At the 
threshold field Eth, the total force acting on a pinned domain of Wigner crystal of linear 
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size £ T is EThPc^{^ T ) 2 - This should be balanced by the total elastic force acting on the same 
domain. Since by definition, the displacement over the length £ T is ~ a, on the order of the 
lattice constant, then the total elastic restoring force must be p m Q 2 a 2 a x a. Therefore, 

p m tt 2 a 2 a 
We can also define j[7| a longitudinal correlation length £ L beyond which the Coulomb inter- 



action energy is smaller than the pinning energy: c L /£ L = c T /(£ T ) 2 , or £ L = 2n(^ T ) 2 /(aa). 
Thus there are two length scales involved in the present problem, related to the shear mod- 
ulus and the long range Coulomb interaction respectively. In the case where the latter is 
screened by free carriers, £ L will be larger than £ T by at least a factor of \j ^fot. 

Various aspects of the FLR model have been studied in the B = charge-density wave 
context [|TTl -|T3|j. The linear response of the pinned state for B / [p,2^] has also been 



determined. In Ref. [|17|], the large difference between the bulk modulus and the shear 
modulus was taken into account. One consequence is the occurence of a /c 1//2 -dispersing 
mode for k lying between 1/£ L and 1/£ T [|5|,|T7||. In the absence of disorder, or for k greater 
than 1/£ T in a statically disordered lattice, there exists a fc 3 ' 2 -dispersing mode which is 
characteristic of a 2D Wigner crystal in a strong magnetic field; it results from the mixing 
due to the magnetic field of the transverse and the longitudinal phonon modes which disperse 
as k and k 1 / 2 respectively when there is no magnetic field, as one can see from the small 
/c-behavior of the dynamic matrix given earlier. The present study focuses on the sliding 
regime, treating the disorder potential as a perturbation relative to the external driving 
force. 

We note that although the present problem bears some resemblance to the sliding motion 
of the vortex lattice in a type-II superconductor ||15| , in that there is a Lorentz force, and 
that the shear modulus is much smaller than the bulk modulus, there is a very important 
difference. The sliding motion of the vortex lattice is driven by the Lorentz force due to the 
superconducting current. The moving vortex lattice then in turn induces a voltage drop in 
the direction of the superconducting current. In the present problem, the depinning electric 



force is supplied externally and the magnetic field enters explicitly the equation of motion. 

The present model may be used to study the sliding of a 2D Wigner crystal with no 

external magnetic field by simply setting the Lorentz force to zero. Such a system has been 



shown to exist for electrons on a helium film |23|| , and more recently, has been argued for 
doped electrons at Si/Si0 2 interfaces |TS[. We will also give our results for cases where 
there is no external magnetic field. When expressed in terms of the shear modulus, they are 



basically the same as those for a sliding vortex lattice [y|. As we mentioned already, the 
difference in the longitudinal phonon spectrum between the two systems does not appear in 
the large velocity limit. 

B. The Hall Resistivity 

Before we embark upon a detailed study of Eq. [I], it is worth noting that one can obtain 
the Hall resistivity p xy immediately from Eq. [I] for a Wigner crystal in the sliding state. 
With a constant current density j x in the ^-direction, the Hall resistivity relates Ey Xt to j x 
via Ey Xt = p yx j x . We disorder-average the Fourier transformed Eq. [I] under the steady state 
condition u x ^ 0. Assuming isotropy, we obtain (using j x = p c u x ): 

p c E e y xt = p m u c u x = pJ — )jx, (4) 

so p yx = B/p c e. This argument is independent of the form of the disorder potential (aside 
from isotropy), and requires only that the current be transported by a k = mode of the 
crystal. We note however that these arguments are strictly only valid at T = 0, where there 
are no thermally activated excitations that carry current. Effects of nonzero temperatures 
on the Wigner crystal state below the sliding threshold are not well understood, and in 
general will depend on the nature of the possible charged thermal excitations in the system. 
The issue of Hall resistivity for the insulating phases in the FQHE regime was raised in 
how to best characterize the nature of the insulating phase |2|| . Kivelson, Zhang, and Lee 
[ 2~3| have suggested that it be better described as a "Hall insulator", which they characterize 
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by a diverging longitudinal resistivity and a constant Hall resistivity at zero temperature. A 
more concrete description of the "Hall insulator" in the FQHE regime has not been given. 
Several authors have recently argued that an Anderson insulator is in fact also a "Hall 
insulator" with Hall resistivity given approximately by B j p c e, p c being the total localized 
carrier density [E1,B3]. 



Experimental results for p xy in the insulating phases have now become available p6 
They are obtained by measuring the differential Hall resistance at a constant current. The 
mixture of p xx into p xy , which can be rather large at low temperatures, can be subtracted 



by an average over p xy (B) and p xy (—B) [pq] . The experimental findings can be summarized 
as follows: 

1). At low temperatures and well below the sliding threshold, the Hall resistivity is given 
by p yx = B/p c , where p c is the total electron density that was introduced into the sample 
during doping, and therefore p xy is insensitive to temperature, in sharp contrast to p xx in the 
same regime which is often thermally activated. When it was measured, the low-frequency 
response was always found to be resistive, with very small out-of-phase component in these 
systems [JSJ. 

2). The Hall resistivity is found to have little temperature or field dependence above and 
below the so-called critical temperature which is marked by the absence of an apparent 
depinning transition. The longitudinal resistivity changes behavior: it may or may not show 
thermally-activated temperature dependence; also it may or may not be ohmic p6| . 



If the applied field is above the sliding threshold for the pinned Wigner crystal, our 
argument above can be used to argue for a normal Hall resistivity. It is not clear at this 
point if the experimental data below the apparent sliding threshold could be understood 
within the Wigner crystal picture. 

It is perhaps worth commenting that at T = the Hall resistivity of a pinned Wigner 
crystal within the linear response regime, i.e., below sliding threshold can also be shown to 
be simply B/p c , following earlier work by Fukuyama and Lee on the magneto-conductivity 



of a pinned CDW |22| . This conclusion can be generalized at a finite frequency below sliding 
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threshold, by repeating the argument given at the beginning of this subsection leading to 
Eq. |], thus avoiding the need to resort to the particular form of the pinned phonon spectrum 



as was done in |22|] . These arguments, as well as those for Anderson insulators |M] , are made 
at T = where there are no thermal excitations in the insulator. However, this mechanism 
necessarily requires that the conduction be carried by the polarization current excited at a 
finite frequency, implying a capacitive response. It then follows that the measured current 
and the voltage must be out of phase by 90°. This scenario is disfavored by the current 



experiments |P6| . More careful measurements geared specifically toward this issue would 
certainly be very useful, however. 

A final comment is in order about the distinction between resistivity and conductivity in 
a disordered system, when the response functions become non-local and local field effects are 
important. In particular, the conductivity tensor is defined by j(q) = J2qi a (<f> ?') ' E{q')\ 
a becomes diagonal only after averaging over disorder. The resistivity is the inverse of the 
conductivity matrix in g-space. The general arguments we have given in this section relate 
to p xy defined as lim^ i^ Q cr~y(q, q') ^ [lim~ -i^VxyiQ, q')} -1 . This is the resistivity measured 
when the current flows entirely in the uniform mode, and may not correspond to the simple 
inversion of a measurement of the conductance. 

C. Steady State Solution with No Disorder 

For VfVdis = 0, the equation of motion Eq. [|, can be solved trivially. In steady state: 

P^E ext = ( " | • v, and , u(f, t) = vt, (5) 

P™ \uj c X J 

with v = Pe E ext . The resistivity tensor p is simply: 

P^f" """). (6) 

and the conductivity tensor a is: 

T _ pI ( x ^ , m 

' Pm (V + " 2 c ){-u; c A ' 
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If there is an additional physical dissipative process which changes A by a small amount 
SX, we will have for the components of the resistivity tensor 

OPxx Pm 



5X Pt 



2 ' 



(8) 



and, 



5p 



•'•// 





sx " Ul 


and for the conductivity tensor 




8&XX 


pI ^ c 2 -a 2 


5X 


p m (cu c 2 + A 2 ) 2 ' 


and, 




Sa xy 


_ p\ -2uj c X 



(9) 



(10) 



5X Pm (cu c 2 + A 2 ) 2 - (11) 

One sees that for A < u c , a xx increases with A, in contrast to the zero-magnetic-field 
case. In addition, \Sa xx \ ^> \Scr xy \ while \a xx \ <C \c xy \ for A <C uj c . 
And finally, 

Sv = -X5X 

v X 2 + uo 2 c y } 

If one views the disorder scattering as an additional source of dissipation which renor- 
malizes A, many of these observations can be carried over. The problem therefore is to find 
how disorder scattering modifies A. To this end, a perturbative analysis is carried out in 
Sec. III. 

III. PERTURBATION THEORY: NONLINEAR CONDUCTIVITY 

We give here the detailed results for the nonlinear conductivity of a sliding Wigner crystal 
in the presence of a strong perpendicular magnetic field with disorder scattering. Results 
here are also used in the next section for studying the correlation functions of the sliding 
Wigner crystal. The dependence of the Hall angle on the sliding current is calculated. In 
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Sec. IIIA, we introduce the response function for the present system, and we outline the 
perturbation theory we use. In Sec. IIIB, we solve for the velocity to the second order in 
disorder scattering, thereby obtaining the nonlinear conductivity. We show in Sec. IIIC that 
the large- velocity perturbation theory also gives a measure of the depinning field. The Hall 
angle is studied in Sec. HID. 

A. The Perturbation Series and the Green's Function 

We shall in this section restrict ourselves to a DC electric field where most of the ex- 
periments are carried out. The extension to the more general DC + AC field case will be 
discussed in Sec. VI. It is convenient in the present problem to write the DC external field 
in terms of a velocity Vq as: 

P^E^=( X ~ Uc )-Vo, (13) 

P™ W A J 

and the displacement as: 

u(r,t) = v t + x(r,t), (14) 

where the first term would be the correct solution to Eq. [I] for the field in Eq. |13] if there 

was no disorder (see Eqs. |6]and[7]). The last term, x, represents the additional displacement 

due to disorder scattering. By substituting Eq. |14| into Eq. |l] and using Eq. [13], we find that 

x(r, t) satisfies the following equation of motion: 

/ A -u c \ , ^ 

p m x + p m \ ■ x + / D(f - r') ■ x(r',t)dr' = - p(f)V '?V dis (r '+ u(r,t)). (15) 

W A J J 

Formally, the solution of the above equation can be written in terms of tensor Green's 

function G as: 

x (r, t)= f dr'dt'G(f- r', t-t')- -p{r')V,pV dis {r' + u(P, t')) , (16) 

where G is the solution of: 
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p m G(r,t) + p„ 



G(r,£) + / D(f-r') • G(r',t)dr' = 6(t)5(r) | ) . (17) 



A -Wr.\ r /l 0' 

UJ C A 

Fourier transforming Eq. [T7|, using the long-wavelength limit expression for D which we 
gave in Sec. IIA, and assuming u <C A, we have: 



G-\k, u)/p m = c L k— + c T k 2 (l - —) - iul 



(18) 



the inverse of which can be taken easily: 

1 



G(k,u) 



P(k,u) 






c L k{l-—)+c T k 2 ^--tu 



k~ 



A u r 



— UJr A 



(19) 



The poles in G are given by the zeroes of P(k, u>) which is: 



P(k,uj) = 27iQ i ap r 



(ka) 3 — k^ — i{ka)K 2 x sgn{uj) 



(20) 



In the above equation, we have defined two dimensionless quantities: 

a; 2 (q; e 2 + A 2 ) 
2-nWa ' 
2 ^A 



(21) 



We have also ignored the a(ka) 2 term whenever it appears together with (ka) term since a is 
small and the relevant (/ca)'s are also small, k^/cl is the wavevector at which the /c 3,/2 -mode 
is at frequency u, and k,\ measures the strength of dissipation at the frequency to. Their 
physical meaning for the sliding state will become more clear later in the paper. 

Eq. [16] can be used to generate a perturbation series for x(f, t) in terms of the disorder 
potential: to first order, x\ is given by Eq. [16] with u(r', t') in Vdis on the right-hand-side 
replaced by vt'; and the second correction x 2 is given by replacing Vdi S {r' + u(r',t'j) by 
xi(?,t) ■ VpVtoi? + vt'); etc. 

B. Nonlinear Conductivity 

Following the last subsection, we obtain for the displacement to first order in the disorder 
potential Vdi S : 
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x(f, t) 



where Fi(k, q, uj) is given by 



d 2 k du d 2 q 

(2^)2 2^ (2tt) 2 



e **-r-**F x (k,q,w) 



(22) 



F x (k,q,u) = -p(k-q) 



G(k,u)-(iq) 



Vdia(q)2n5(u + q- v); 



(23) 



and for the velocity to second order in Vdis'- 

A -u c 



Pr, 



where F 2 (k, G) is: 



U0 r A 



<*> = £Ip(g)I s 



G 



d 2 k 



F 2 (k,G), 



(24) 



F 2 (k) = kT{k) 



k-iG(-k + G,k-v) ■ k 



(25) 



In Eq. [24], we have averaged over disorder and integrated over the 5-function. Vdi S {q) = 
J d 2 re~ l ^ r V ( n s (r), and T(k) is the disorder correlation function: < Vd is (k)Vd is (q) >= 6(k + 
q)T(k), where < ... > indicates averaging over disorder. We have also used the fact that for 
an undistorted Wigner cyrstal: 



p(q) = (2ir) 2 J2p(G)5(q-G), 
G 



(26) 



where G's are reciprocal lattice vectors. 



Using Eqs. |23] and [^ in Eqs. |22] and [24], we have: 



and, 



Pr. 



x(r, t) 



A -a> n 



LU r A 



r+i(q-v)t 



(2tt) 2 (2tt) 



-p(fc - g) 



G(k,-q-v) ■ (iq) 



V dis (q), 



(27) 



(•*' 



} = Elp(^)P 



G 



d 2 q 
(2tt)' 



qT(q) 



q- (— JmG(— g + G,q-v))-q 



(28) 



In both Eq. [27j and Eq. ^, the intergration over k is over the first Brillouin zone, and 
that over q is over the whole reciprocal lattice space. Eq. ^7J is used later for studying 
various correlation functions in a sliding Wigner crystal. We shall now focus on Eq. B3 
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which gives the lowest order correction to the conductivity. It is most easily studied using 
the representation of G in Eq. [19|. 

Before we proceed further, we should mention that an alternative view of the same 
perturbation theory is to ask: In order to sustain the current flow p c v in the presence of 
Vdisi what additional external field one would need above that given by Eq. [H| The answer 
is that to second order in disorder potential the additional field E p is given by: 

Pc EP = J2\ P (G)\ 2 Jj^L qT(q)\q.I mG (-q + G,q-v)-q\. (29) 

The integral in Eq. [28] or Eq. ^ is dominated by two infared terms where the elastic 
medium theory is expected to be valid. They are respectively related to ^ and k\ in Eq. |20| . 
For the convenience of the reader, we rewrite them here at u = u v = \G ■ v\: 

c^ + A 2 ) 
ft,., — 



2vrfi 4 a ' 

We have suppressed the dependence of u v , n^ and k\ on G for ease of writing. Both ^ and 
k>\ must be much less than unity in order for the elastic medium theory to be valid, and 
their ratio k w /k\ ~ (cj 1 ,o; 4 a/47r 2 A 3 f2 2 ) 1//6 depends on the sliding velocity v through u> v . The 
dependence is not very strong, ~ v 1 ' 6 . 

The result of the integral in Eq. ^B| depends on the relative size of n^ versus K\. We find 
that in general: 

C I • (£) = -a — J2 A(G)G sgn(G ■ v) 
u c A / Pm G 

= -ocdu-, (31) 

V 

where 

A(G) = -^\p(G)\ 2 T(G)\G\\ (32) 

The constant «o in Eq. |3I]is tt/3 if k^ ^> k\ and is 7r/4 if k w <C k\, and changes monotonically 
between these limits. a<&s is uncertain up to a numerical prefactor of order unity that 
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depends on the direction of v relative to the lattice orientation. The relative stability of 



various sliding directions was discussed in [Rj in the sliding vortex lattice context from a 
minimum-entropy-generation argument. In this work, we simply assume that the Wigner 
crystal is sliding along one of the stable configurations which may be defined with a vector 
G. 

Therefore, if we use the alternative approach given by Eq. |29], we find that for a given 
average sliding velocity v, the total applied field E ext must be: 

A —OJr. 



E 



ext 



Pn_ 
Pc 



Ur 



■ v + a dis - 
A l v 



(33) 



It is important to note that the right- hand- side of Eq. |3l] does not depend on the mag- 
nitude of v, only on its direction. This feature is common to all the 2D systems studied 
within the second order Born approximation, and has been noted in the CDW context in, 
e.g., Ref. |12|. Eq. [U] also shows that the net effect of disorder is to increase the A in the 
equation of motion by an amount of a^, s jv . 

It is straightforward to obtain the total velocity v + (x), for a given external field E ext 
in Eq. O. The resultant conductivity tensor has components given by: 



a , 



PI: 



p m A 2 + u\ 



A + 



Ctdis U c 



A 2 



v w 2 + A 2 



(34) 



and 



a 



■>-!l 



— O ni 



Pc 



iO r 



ctdis 2u; c A 



yx ' n \2 _i_ , ,2 

Pm A + ^ c 

Or equivalently in terms of the resistivity to order a^ g : 

\p 



v cj 2 + A 2 



(35) 



Pi 



Pi 



1 + ctdis/ {vX) 



(36) 



and 



Pxy 



UcPr, 
Pi 



B/p c 



(37) 



Notice that p xy is not modified to order ctdis, as was expected from the general analysis given 
earlier. The above equations are the same as Eqs. |q to |lT| if 8X — adi S /v. 



The v appearing in Eq. |3T] is the magnitude of the average sliding velocity. If one does 
not require self-consistency in the solution for the sliding velocity in the presence of disorder 
potential, one then must take v = vo from Eq. [13]. However, Eq. |31| can also be solved 
self-consistently for the velocity, giving a self-consistent v. These two procedures give the 
same result to order a^s, but differ in the regime where the external field is not much larger 
than the sliding threshold. The alternative approach given by Eq. |29| starts off with the 
self-consistent velocity, and thus automatically gives results that are self-consistent, and is 
in fact identical to what one obtains by solving Eq. |3l] self-consistently. However at fields 
not too far from the threshold field, a priori, it is not clear which one gives the better 
solution. We shall come back to this when we discuss the Hall angle, and give more details 
for results at fields not too larger than the threshold field where these two treatments differ. 
Our results for the conductivity and resistivity above were given to order a^s, at which level 
the self-consistent theory and the non-self-consistent theory agree. 

Our results for the conductivity and resistivity can be used to compare with experiments. 
In Fig. 1, we show the calculated I — V characteristics of a sliding Wigner crystal in the 
presence of disorder scattering. We assume a Hall bar geometry of measurement, i.e., we 
assume a fixed /, and plot the longitudinal (Vl) and Hall (Vjr) voltages as a function of 
/. In Fig. 2, the differential resistivity is depicted as a function of the applied voltage, also 
in the Hall bar geometry. Our results are not expected to be reliable close to depinning 
threshold. In particular, the exponent of the I — V curve given by the present perturbation 
theory is incorrect near threshold. The dotted portions of the lines are meant to convey 
the approximate nature of our results in this regime. We shall compare these with the 
experimental findings 0-0] later in this paper. 

C. Relation to the Static Depinning Field 

We here show that the quantity p m o>disl Pc appeared in the last subsection is also a 
measure of the threshold field which has been defined in Eq. 2 in terms of properties of the 
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static Wigner crystal. The static Eth, from Eq. 2, is fL7j: 



p m c T a 



2 ' 



fr~*=^, (38) 



where n, and V =< V 2 >*/ 2 are the average impurity density and strength. Combining 
Eqs. [SB] and putting c T = aQ 2 a 2 , we have 

riiV 2 /p c a 3 , . 

E Th ~ o22 - 39 

p m mz^cr 
We may also estimate Exh as the field at which the second-order correction to velocity 



becomes comparable to the zeroth-order term. From Eq. 31, we estimate 



E ^ adispm = y- \p(G)\ 2 T(G)/ Pc a 3 
pc q p m aQ 2 a 2 

Indeed the above two expressions are equivalent. Therefore the high- velocity perturbation 
theory in fact gives an estimate for the threshold field which can be used to estimate other 



quantities such as the correlation length as was done in Ref. |fL7|| . Using this, one may also 
relate the quantity A(G) defined earlier to the depinning threshold field. We then have: 

^a J2\G\A(\G\) = Pc E Th . (41) 

G 

D. Dependence of Hall Angle on Sliding Current 

As we mentioned already, the relationship between the external field and the total sliding 
current depends on whether one requires self-consistency in the velocity for fields not too 
large compared to the depinning field. 



In the case where we do require self-consistency in Eq. 31, i.e., the velocity that we use 



on its right-hand-side is not v in Eq. O, but the self-consistent v, we find: 



-a dls \ + y^LA 2 + (u; 2 + X 2 )(p 2 Ey P l - a 2 dls ) 

^ + A 2 • l42j 
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For E 3> ocdisPm/Pcj it reduces to p c E/(p m Juj 2 + A 2 ); and for E only slightly greater than 
OidisPm/Pc it is (p c E/p m — a.dis)/\. The Hall angle at an arbitrary velocity v is: 



tanO 



LO c V 



H 



(43) 



Xv + pmOidis/Pc 

It becomes 45°, i.e., the field along the current becomes equal to the field perpendicular to 



it, when: 



A + a d i S /v = uj c 



(44) 



Putting in the velocity in Eq. |^2|, substituting a dis with p c E Th /p m as we argued in the 
previous subsection, and solving for E, we find: 



-£/45° nz UJ r , 



V2- 



(45) 



Eth oj c — X 

For E less than that given in Eq. ^, the Hall angle is smaller than 45° and the field and 
the current are largely parallel. It becomes greater than 45° for greater E. For uj c less than 
A, the Hall angle can never reach 45°. 

The results are different if one does not require self-consistency in the velocity, i.e., if 
one uses v from Eq. [13] in the right-hand-side of Eq. ^I[ In this case, the velocity is related 
to the external field as: 



\ 



E 2 + E 2 Th - 2EE Th X/JX 2 + u, 



(46) 



A 2 + u 2 c 

where we have used Exh for p m adi S / p c - The large velocity limit remains the same, but in 
the limit of E only slightly greater than Eth, it becomes, 



3 



Th* 



2(l-A/ V / A 2 +o; c 2 ) 
\l A 2 + co 2 



(47) 



Now the Hall angle is: 



tanO 



uo c E - 2Xuo c E Th /JX 2 + u t 



H 



XE-(X 2 -u 2 )E Th /JX 



(48) 



ui. 
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It becomes 45° when the total velocity is: 



^cErh 



v & = V2 ^ T t h (49) 

We should add that these results must be treated with caution at fields just above the 
threshold. Both approaches given above are extrapolations of the perturbation theory to 
a regime where perturbative corrections are no longer small compared to the unperturbed 
term. However, if we identify E Th = a,n s p m / p ci then the self-consistent solution is more 
sensible for applied fields not much greater than Etk = O-disPml 'p c , in that it gives a sliding 
velocity v oc {E — ETh)- Although the exponent relating the velocity to the field is incorrect, 
v does approach zero at the depinning threshold field Eth- 

IV. CORRELATION FUNCTIONS AND LOCAL STRAIN 

To better understand the state of a sliding Wigner crystal, we consider in this section 
some correlation functions. In charge- density waves, examination of the strain has proven 



useful in understanding the qualitative features of the Fukuyama- Lee- Rice model \\ I] and 
determining the limit of validity of the elastic model. We use the perturbation theoretical 
analysis in Sec. Ill to study the correlation functions of the magnetically-induced Wigner 
crystal under the sliding condition. We show our results for the velocity-velocity correlation 
function in Sec. IVA. A more detailed study of the correlation length in the sliding state is 
given in Sec. IVB. The results for the strain are given in Sec. IVC. 

A. Correlation Functions 

In the sliding state, only the fluctuating parts of the displacement and velocity need to 
be considered in the correlation functions. We begin with the solution to Eq. [15| for the 
displacement in the DC limit: 

^' t) = I ^ e ^ (_p( ^ ~ ® [ G( ^' ~ g "' ^ ' N VdiM - (50) 
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x(k, t) is the Fourier transformation of x(f, t). 

After averaging over disorder, one has for k <C G for equal-time correlation function: 



< x(k,t)x(-k,t) >=Y,\P(G)\ 2 T(G) 

G 

and similarly for the velocity, 



G(k, -G-v)-G 



G(-k,G-v)-G 



(51) 



< v(k,t)v(-k,t) >=Y,\G- ^| 2 |p(G)| 2 r(G) 



G 



G{k, -G-v)-G 



G(-k,G-v)-G 



(52) 



In perturbation theory on the FLR model, Matsukawa and Takayama [13|] and Fisher |2T 



have argued that the expression for the positional correlation function contains divergent 
terms at 4-th order in the perturbation series, but the expression for the velocity correlation 
function does not up to the highest order terms examined. For this reason, we view our 
velocity-velocity correlation function as being more reliable. However, at the present level 
of perturbation theory, the two correlation functions behave in almost exactly the same way 
and are described by the same correlation length that we discuss in the next subsection. 

The correlations at large distance can be obtained by the Fourier transformation of the 
above equations. The r-dependence of the real-space correlation functions is given by the 
/c-dependent part of the fc-space correlation functions: 

d 2 k 



g(r) 



{%*)< 



e lk - r \G(k, -G ■ v) 



(53) 



The imaginary part of the poles in the Green's function G determines the characteristic 
lengths in these correlation functions. For these dynamic quantities, the magnetic field nec- 
essarily enters in contrast to the static case. We discuss these lengths in the next subsection. 

B. Dynamic Correlation Length 

Before we describe our results for the correlation length which is contained in Eq. |53|, 
it is useful to review the results for a CDW system treated within the FLR model. The 
corresponding equation of motion was given in Eq. |2|. In the CDW case, the Green's function 
takes a simpler form: 
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G CD w{k,uj) = — - ; — — . (54) 

Ak 2 — ip m AUJ 

Here A is the CDW elastic constant. 

The correlation functions for the FLR model in 3D have been calculated by Matsukawa 
and Takayama ||13| . One can easily extend their result to ID and 2D. We find that the 



velocity-velocity correlation as defined in the last subsection in the large r-limit becomes 

r 

{-h=) D ~ l e iv ^ times a purely oscillatory factor sin(— 77= — 0o), where 0o is a phase factor 
which depends on dimensionality D, and the correlation length £„ is given by 



& = y/A/( Pm XG\v z \). (55) 

This expression for the correlation length holds as long as u <C A, i.e., as long as the 
co> 2 -term can be ignored in the Green's function so that the oscillation at uj is purely diffusive. 
£u is nothing but the decay length of an oscillation at frequency uo = \G ■ v\. It describes the 
distance over which a disturbance at frequency ui v will propagate. 

This "dynamic correlation length" can be related to the static Fukuyama-Lee-Rice length 
£ in the CDW case as: 



V ¥p m \G\v z \ \jE-E Th 

So for fields E not too large compared to the sliding threshold field E^h, the dynamic 
correlation length is comparable to the static correlation length in the B = CDW case. 

With this picture in mind, we examine the present situation which is more complex. Since 
it is clear from the above arguments that the dynamic correlation length is determined by 
the excitations of a Wigner crystal at the washboard frequency u v , we examine the excitation 
spectra which appear as poles in the Green's function, given by Eq. |20] and Eq. ^. Setting 
the denominator of the Green's function to zero, one has simply in terms of k^ and k,\: 

(kaf -Kl + i(ka)K 2 x = 0. (57) 

We are interested in the case where k u and k\ are given by those at uj = u v = \G ■ v\. 
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Note that Eq. [57] is a cubic equation for k. So for each given u, one finds generally three 
poles in the complex fc-plane. If u> v is very small, then the modes are once again diffusive. In 
this frequency range, there is no propagating g 3 / 2 - dispersive mode. More quantitatively, at 
very small average sliding velocity, for k w < kai the three roots of Eq. |57| are approximately: 

k\ 2 a 2 + %k\ = 0, (58) 

and, 

- kI + ik 3 aK 2 x = 0. (59) 

These poles give rise to two lengths, as: 

6,2 ~ -j — ~ — , (60) 

"4,2 K X 

and, 



K 



6 ~ h ~ a i' (61) 

The largest one of the three lengths, £3 , will determine the large- distance behavior of the 
velocity-velocity correlation function, although stronger correlations between the velocities 
of different elements will set in when their distance becomes comparable to £1^. The first 
one of the two lengths (£1,2) in this over-damped limit is similar to the CDW case, since: 



— = yl<?/{ Pm \\G-v\). (62) 

K x 

In the limit of E only slightly greater than Efh-, it can be rewritten using the small w-limit 

of Eq. H as: 



- = W-ET^- (63) 

k x \E-E Th K J 

But in the regime where E is much larger than Etk while k^ is still much less than kaj we 
have instead: 



i^T^ (64) 
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which differs from the CDW result in Eq. [5(| 

As the sliding velocity increases, we get into the regime where n u ^> K\, The three roots 
are now approximately: 

K 2 

k x a = k^ + ir-^-, (65) 

and, 

k 2t3 a = ^e ±i27r/ \ (66) 

In this case, the /ci-mode is almost propagating. The corresponding lengths are related to 
the static correlation lengths as: 

6 ~ a-j, (67) 

K \ 

and, 

6,3 ~ — • (68) 

In this limit, £i still describes the decay length of an oscillation at frequency u v , but the 
other two describe the dispersive propagation within a length that is smaller than the decay 
length. This is due to the particular dispersion of the /c 3 / 2 -mode. This phenomenon does not 
appear in the usual CDW systems where excitations at u v are always overdamped as long 
as lo v <C A. In the present case, the magnetic field makes it possible that even for w„< A 
in the sliding state, there is still an almost propagating mode. It also has interesting effects 
on the strain in the sliding Wigner crystal, which is the subject of the next subsection. 

For a discussion of the crossover between the two regimes in a typical experimental 
system, we refer the reader to Sec. VIIIB. Under most experimental conditions where the 
DC current is a few nA, k u and K\ are comparable. Since their ratio has only a weak 
dependence on the sliding velocity (~ v 1 ^ 6 ) for all cases of practical interest, neither limit 
that we discussed above for the velocity-velocity correlation length is strictly applicable. 
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C. Strain 



One can also calculate from Eq. 50 the local longitudinal (L) and transverse (T) strain 



in the sliding Wigner crystal, defined as: 

E s L =<(V-x(r,t)) 2 >, (69) 

and, 

^=<(Vxf(f,i)f>. (70) 

Eq. [F(] gives the local elastic energy associated with the shear modulus; and Eq. |69] would 
be the local elastic energy associated with the bulk modulus, if there were no long-range 
Coulomb interaction. 

From Eq. [H], we find that: 

El = J2\p(G)\ 2 T(G) j ^ 2 \k-G(k,-G-v)-G\\ (71) 



and, 

Et = Y, \P(G)\ 2 T(G) J-^fx (G(*. -G ■ v) ■ G)\ 2 . (72) 

Ej, can be transformed into 

E> T = Z 1,(0)1^(0)10^ ^/; d(M - W' (73) 



G 

Once again, the result of this integral depends on the relative size of k^ and k x , or, more 
physically, depends on if there is an almost propagating mode at the washboard frequency. 
If not, i.e., for k^ <C K\, 

E strain = £ \ p{ G)\ 2 T{G)G 2 ± In- = P^-ln-, (74) 

^ 47rp4(c i ) 2 k x il 2 aa k x 

while for k^ ^> k,\, 

E stroin = j2 \ p{ G)\ 2 Y{G)G 2 — J— -{In— + |4) = P7^( ln - + 1% ( 75 ) 
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Here (3 is a dimensionless constant of order unity, and aidis is ~ p c ETh/p m , as was given in 
Eq. |31]. As we have mentioned already and will discuss further in Sec. VIIIB, the second 
term on the right-hand-side in Eq. [75] is smaller than the first in the experimental systems. 

The first term in Eqs. |73| and |75] diverges as v — > as In ~ Th . The second term in 
Eq. [75], which is due to the fci-mode in Eq. ^5], has no counterpart in the usual charge-density 
wave cases. This is due to the fact that in the CDW cases, one only considers the strongly 
diffusive mode as long as the washboard frequency is smaller than A. In the present case, 
however, under the same condition, there is the possibility of having a mode like k\ in Eq. |65| 
whose real part is much larger than its imaginary part. It directly reflects the fact that the 
/ci-mode is only weakly damped, causing an almost resonant absorption at (oo v , ki). 

Within our high velocity perturbation theory, the local transverse strain divergees weakly 
as v — > 0, suggesting that the elastic medium theory used here will break down before the 
pinning transition is reached. However, it is possible that higher order terms may change 
the form of the divergence found in Eq. |7^ and Eq. [75], or eliminate it entirely. In all cases 
^stram j g sma ll e r than E^ ram by at least a factor of (27r/a) 2 , and does not contain the 
divergences that we discussed above if the Coulomb force is unscreened. 

The possibility that the elastic model breaks down at small velocity has been discussed 
by Coppersmith and Coppersmith and Millis fl4| for the Fukuyama-Lee-Rice model without 



a magnetic field. Using Eq. [74] or Eq. |75] and typical material parameters for a FQHE device 
outlined later in this paper, we find that the transverse strain at a typical transport current 
of InA is slightly less than 10%; the longitudinal strain is completely negligible. 

The static correlation length is determined by the shear modulus alone, because the 



shear modulus is much smaller than the bulk modulus |fL7| . We see that in the sliding state, 
disorder scattering also induces principally transverse distortions. This is confirmed in both 
the study of the nonlinear conductivity in the previous section, and here by examination of 
the strain distribution. However for the velocity-velocity correlation length, the magnetic 
field and the bulk modulus both enter. 
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V. FREE CARRIER EFFECTS: LARGE VELOCITY LIMIT 

In the previous sections, we have studied some aspects of the sliding motion of a Wigner 
crystal in the fractional quantum Hall regime. We now consider the effects of free carriers, 
which may be characterized by a separate conductivity tensor, on the sliding motion of the 
Wigner crystal. These free carriers may be thermally excited, or due to sample inhomo- 
geneities. We treat them phenomenologically. 

We assume that the effects of the free carriers can be represented by a conductivity tensor 
whose diagonal and off-diagonal elements are (<t xx , ±cr xy ). Under the assumptions that 

1) there is no interconversion between the free carriers and the Wigner crystal, and, 

2) the coupling between the free carriers and the Wigner crystal is purely electrostatic, 
it can be shown [17 that longitudinal component c L in the Green's function becomes: 



c L k -> fl 2 



+ (l + a)(ka) 2 



1 + ika^ 



tjQUJ 

o2 



(76) 



Here ero = ae£l/2ir, slightly smaller than e /h for a typical GaAs/AlGaAs modulation doped 
FQHE sample. 

It is therefore clear that the importance of free carriers is measured by the magnitude of 
the quantity Mp = kaa xx Q/aoUJ. We are interested in the range of (k, ui) which corresponds 
to that of the disorder scattering potential. The characteristic frequency of the disorder 
potential is u v ~ v/a in the sliding state of sliding velocity v. For a current of a few nA, 
which is typical of most experiments, v is ~ 1CT 2 m/s. This yields uo v ~ 10 6 Hz, much smaller 
than Q. This has two consequences. The first is that for such frequencies G(k,u>) can be 
approximated by that given by the elastic medium theory, which is true with or without 
the free carriers. Secondly, combining this with the phonon dispersion which determines the 
characteristic value of k, we conclude that Mp ^> 1 as long as o xx Joq of the free carriers is 
greater than 10~ 3 , as a conservative estimate |27[ . 



This range of a xx in fact covers most of the experimentally measurable range. We 
therefore conclude that for our purposes, we only need to consider the case Mp 3> 1. We 
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then have: 



c L k = -iwAi + Vt\ka)\ 



with 



Ai = 2vrfi 



CTq 



<Ji 



and the corresponding Green's function is given by: 
G-\k, u)/p rn = n 2 a 2 k 2 ^ + aQ 2 a 2 k 2 (l - ^) 

rv rC 



, kk 

— iuAi— — iuo 



A -u c 

UJn A 



Its determinant P(k,u) is 



P(k,uj) = aQ A (kaf-u 



uj 2 + A(Ai + A) 



2/t>„\2 



iuVt (ka 



a\\ + (1 + a)A 



(77) 



(78) 



(79) 



(80) 



From here most of the results from the previous sections follow without major changes. We 
therefore only list the results in the rest of this section. 

A. Nonlinear Conductivity 



The perturbative calculation proceeds as in Sec. III. We find it convenient to define in 
the present case three dimensionless quantities: 



/v, 



\LU V \U) C 



(81) 



\LO 



n 2 



m+A/a), 



(82) 



Q 2 



(Ai + 2A), 



(83) 



where u v = G ■ v is the sliding washboard frequency, and k% is less than k 2 . K\ja is the 
wavenumber at which the lower hybrid phonon mode, which disperses like uj 2 = aQ 4 (ka) 4 /u 2 
in the present case, matches the washboard frequency, and k 2 and % measure the strength 
of dissipation in the present context. 
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In the limit of uj c ^> A and Ai 3> A/a, we find that the spatially averaged velocity 
correction to second order in disorder scattering v 2 is: 

r) -^ = ~ i to QW Y.\p(G)\ 2 nG)\G\*G sgn{G-v) = -«* (84) 

This yields a resistivity tensor which is identical to that in the absence of free carriers. 
This is because the leading corrections depend solely on transverse distortions, and only the 
longitudinal modes are affected by free carrier screening. The results were given in Sec. Ill, 
and will not be repeated here. At higher orders of the perturbation theory, the longitudinal 
channel affects the results; and the conductivities with or without the screening effects will 
be different. 

B. Correlation Functions 

As in Sec. IVB, the correlation length of the sliding state for all of the correlation 
functions will be determined by the poles of the Green's function G(k,u), given by: 

(kaf-K\ + i(ka) 2 Kl = 0. (85) 

For simplicity, we assume k±^$> K2, and obtain: 



e = 4a/t 1 /^~l/ v /|^|, (86) 

which approaches infinity as \oo v \ — > 0. 

The evaluation for the strain also proceeds as before. We here only give the results: 

E s H v x , T = E W«i/(Wg«, (87) 



and 



ei = ( v . *)> = e \*w&J^J(*mw {( Jw+°* (88) 

We see once again that E§, ^> Z?£, although the actual ratio depends on the details of 
the relevant parameters. For K\ ^> K2, we have: 
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E§ = £ |p(g)| 2 r(G) , '^f z 9 Zn- = ^Zn-. (89) 

Here f3 is once again a constant of order unity. The logarithmic divergence has the same 
origin as that for the transverse strain distribution function previously. It is the result of 
the dimensionality of the present problem. 

VI. AC + DC INTERFERENCE EFFECTS 

We now consider the possible interference effects due to the presence of both a DC field 
and an AC field. It is more convenient here to adopt the alternative approach mentioned in 
Sec. IIIB and determine the electric field E Tot required to sustain a current density j which 
has both DC and AC components, given by: 

j = PcVDC + p c v A ccosuj Ac t, (90) 

within the perturbation theory, 

E lot = e + Ep, (91) 

where 



A —u, 



c 



p c E = p m \ ) (VDC + VACCOSUJAct), (92) 

L)r A 



and to second order in disorder potential 

d 2 q 



p c E p = -i J2 \P(G)\ 2 j ——^qr(q)J ni (q- v A c/u A c)Jn 2 (q- vac/uac) 

<3,ni,n 2 ^ ' 

e iu A a(n 2 -n l)t /-,. Q( _^ + £ ^^^ + -.. J DC ) . f\ . (93) 



The notations here are the same as before. J n (x) is the n-th order integer Bessel function. 
The sums over m, n-i include both positive and negative integers. The derivation of Eq. S3 is 



the same as its counterpart without the AC field and essentially follows that in e.g. Ref. JT2 
Therefore it will not be repeated here. 
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There are two inteference effects: 

1). For rii = n 2 , the resultant DC field is modified by the AC current (Shapiro anomaly); 

2). The linear AC-response, i.e., the cuAc-component of E p , will be affected by the DC 
current. 

Both effects enter through the frequency dependence of the Green's function G where 
nujAc + Q • vdc appears. We discuss them separately in the next two subsections. 

A. Shapiro Anomaly 



Setting ri\ = n 2 in Eq. [53], we have: 

d 2 q 



/nn / - \ 

-—-qT(q)(J n (q- v AC /u AC )f I q ■ G(-q + G, nu AC + q- v DC ) ■ q) ■ 

(94) 



(2tt) 5 

G,n 



The n = term in the above equation gives us the nonlinear conductivity in the absence 
of the AC field, which was studied in previous sections. If n ^ 0, under the condition that 
\G ■ vdc + nu)Ac\/\G ■ vdc\ 3> max(K\, kJ), i.e., we are somewhat away from the resonance, 
we find: 

Pc & = a J^ Gsgn{G ■ v DC - nu AC )J 2 n {G ■ v AC /u AC )A{G). (95) 

G,n 

A(G) here is the same as that given in Eq. [32[ 

We here fix u A c to be positive, and have used the fact that in summing over G, ±G 
enter together. a in Eq. |95| is once again 7r/3 if n u ^> k x and 7r/4 if k u <C k a . 

If we further relate A(G) to the threshold field Etk through a^ is as was done earlier 



in the paper in Eq. |32| and Eq. |4T], we finally have for a DC current j = p c voc, the total 
external DC field is: 

/ A -u c \ 

p c E = p m \ • v + PcEtiiVdc/vdc 

\u c X J 

+ X ^ \G\A(G)J n {G ■ v AC /u AC )sgn(\G ■ v DC \ - nu AC )v DC /v DC . (96) 

G,n 
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In cases when G is not parallel to v^ci as long as v^c is along one of the six nearest neighbor 
or the six next nearest neighbor G- vectors, summing over all G"s with the same G-voc would 
yield a correction to the field in the direction of vdc- 

We therefore see that the Hall resistivity is once again not altered. The longitudinal 
resistivity shows upward jumps with increasing current, or downward jumps with increasing 
frequency of the AC field. If one is instead interested in the conductivity, the direction 
of the jumps will depend on whether or not the Hall angle is greater than 45°. If not, 
downward jumps in resistivity will translate into upward jumps in conductivity; otherwise, 
downward jumps in resistivity will become a downward jump in conductivity. This is familiar 
in experiments in the FQHE regime |I] where one often encounters a situation in which the 
off-diagonal elements of the conductivity /resistivity tensor are greater than the diagonal 
elements. The maximum value of the G-vectors that can give rise to the jump is cutoff 
by either the G-dependence of disorder correlation function T(G) or by the density p{G) ~ 
e -\ G \ 1 b/ 2 . Both enter through the expression for A(G). 

This is the Shapiro step in the DC voltage, on the order of the depinning field but 
smaller in magnitude, appearing when we vary either the frequency of the AC field or the 
DC current. (p c E p is proportional to the DC voltage.) The main experimental signatures 
are therefore: 

1). If one sweeps the frequency while keeping the DC current constant, as uuac crosses 
\G ■ vdc\ from below, the DC voltage decreases, leading to a decrease in the DC resistivity. 
The jump in DC voltage should be proportional to the threshold field, with the proportion- 
ality coefficient being J^{G ■ vac/uac), which gives its dependence on the amplitude and 
the frequency (or the DC current) of the AC field. 

2). If one instead gradually increases the DC current at a fixed AC frequency, one 
observes an increase in resistivity. 

3). In usual CDW systems, one often measures the differential resistivity and looks for 
peaks at resonance. The width of the peak gives additional information on the correlation 



34 



length. In the presence case, the situation is more complicated and the width will depend 
on the velocity-velocity correlation length which we discussed extensively earlier. 

Overall, neither p(G) nor T(G) decays very fast with G for small filling factors with small 
magnetic length. There will be many more peaks than in the traditional CDW systems where 
p{G) is finite at a single G and its inversion. The features, if indeed observable, are expected 
to be much broader because the contributions come from many G's. Perhaps a more severe 
hindrance is that the correlation length in the Wigner crystals around v = 1/5 appears short 
17| . Traditionally, it has been easier to observe the interference effect in the AC response. 



We turn to it in the next subsection. 

B. Linear AC Response Function 

To obtain the linear AC response in the presence of a DC current, we assume small vac 
and only keep terms linear in vac i n Eq. |93[ We have: 

PcEP = -j^y 2 e-— t j:\p(G)\ 2 

J d 2 qqT(q)(q-v A c/2u AC )(q- (-G(-<f + G, q ■ v DC ) + G(-q + G,uj A c + <f ■ v DC ) • q) 
+c.c, (97) 

from which we can determine the in-phase and the out-of-phase component of the AC resis- 
tivity. For the in-phase part that is proportional to cosuiAct, we find: 

(p c EP) in = a ^2 Gsgn(u AC ~ G ■ v DC )(G ■ v ac /uac)A(G). (98) 

G 

And for the out of phase part that is proportional to sinuiAct, we have: 

(Pc&) out = £ GA(G)(G ■ v A c/uJAc)ln max{ ^'' Kx 'l (99) 

where k w and K\ are defined as before and the primed quantities are obtained by replacing 
G ■ vdc with ujac — G ■ vdc in the respective unprimed quantities. We have in terms of the 
threshold field: 
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( Pc EP) m = Pc E Th £ G/\G\sgn(uj AC - G ■ v DC )(G ■ v AC /u AC ), (100) 

G 

and 

{p c EP)out = Pc E Th /4a J2 G/\G\(G ■ v AC /u AC )ln maX ^'' Kx '} , (101) 

where the sums are over those that give the same G ■ v among the six smallest G's only. The 
above equations also hold for larger G's, with the necessary replacement of Etk i n Eq. |100| 
by an expression in which the |G|-dependence of various terms in A(G) is restored. 

The out-of-phase part of the AC resistivity should show an inductive to capacitive tran- 
sition as u A c goes from below ~ G ■ v^c to above. For large ujac, it becomes -^—lnuj A c- 
We also note that it has only a component along the direction of Vpc, that is, it only enters 
the longitudinal resistivity. 

The in-phase term involves now three vectors: G, the DC current vdc, and the AC 
current v A c- In cases when the latter two are not aligned, it will enter both the diagonal 
and off-diagonal AC resistivity. In particular, if the relevant G satisfies G ■ v A c = 0, the 
jump in the real part of the AC resistivity disappears, although the imaginary part will still 
show an inductive anomaly. If v A c is parallel to vbc-, then 

(p c EP) in = p c E Th sgn{u AC - \G ■ v DC \)/{au) AC )v Ac . (102) 

The equation is uncertain up to a numerical factor smaller than unity. The diagonal re- 
sistivity changes from roughly %(A — adi S /au A c) to ^r(A + otdi S /au} A c) as uj A q increases 

Pc Pc 

through G ■ vdc- Approximately, for the out-of-phase part, we find: 

(p c E) out = )^A(G) In = . (103) 

^ o-^ac G ■ vdc 

We also note that when the DC current vnc an d the AC current v A c are not parallel, 
the system no longer exhibits the rotational invariance that we invoked in arguing for a 
classical Hall resistivity. Indeed, in this case, we find that the Hall component of the in- 
phase AC resistivity will not be B/p c , with corrections occurring first in second order in 
disorder potential. 
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Care must be exercised when converting these statements about the changes in the AC 
resistivity to those in conductivity, as in the case of the Shapino anomaly that we discussed 
in the last subsection. When the Hall angle for the AC current is smaller than 45°, upward 
jumps in the real part of the diagonal resistivity translates into downward jumps in the real 
part of the diagonal conductivity, and the inductive to capacitive transition remains. But 
for Hall angles greater than 45°, the situation will be reversed, we therefore choose to give 
us results in terms of resistivity. 

VII. RESULTS FOR ZERO MAGNETIC FIELD 

In the case of a zero magnetic field, the Green's function has the same form, with u c now 
equal to 0, and P{k,uj) replaced by: 

P(k, u) = 2naVL A p m [(kaf - K 3 W - iK, 2 x {ka)sgn{uj)], (104) 

where as before, we have ignored the a(ka) 2 /a7r term in the last term on the right-hand-side, 
and 

4 = j^, (105) 

and, 

All of the results for the resistivity are in exactly the same form as those of the diagonal 
resistivity with the magnetic field, if we set u c = and put in the new k\ and k u . We will 
not repeat them here. The only exceptions are the velocity-velocity correlation length. We 
give the new results here: 

1). For k w <C K\, i.e., the overdamped case, 



«u = « T ^H & = £VV ; (107) 

V hi — &Th h, — £j Th 



2). For Ku, > Kx, 
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6,3 = e$) ll \^-?'\ 6 = ^{%f\^-) 1,z . (108) 

4 & — bj T h 4 & — bj T h 

In the overdamped region as specified in case 1, the equations are the same as for the 
CDWs. Case 2 is where the differences lie. 



VIII. SUMMARY OF OUR MAIN RESULTS 

In Sec. VIIIA, we summarize the main theoretical results regarding the transport prop- 
erties of a sliding Wigner crystal obtained in this work. In Sec. VIIIB, we estimate various 
characteristic lengths and frequencies that have appeared in our theory using parameters 
appropriate for a typical high mobility modulation doped GaAs/AlGaAs heterojunction. 
Our calculations should be strictly valid only in the very low temperature limit, while ex- 
periments are done at a variety of temperatures. Therefore in Sec. VIIIA, we discuss the 
possible temperature effects on the nonlinear I — V curves. 

A. Summary of the Theoretical Transport Properties 

We summarize our main results in the form of predicted I — V characteristics of a sliding 
Wigner crystal. 

A). Hall bar geometry 

In the Hall bar geometry, one controls the current that runs through the sample. If 
j = p c v along the x-direction, then, 

E x = ^j x + E Th , (109) 

Pc 

Ey = ^J X . (110) 

Pc 

We notice that the functional form of Eq. |109j is in very good agreement with the ex- 
periments by Li et al. || and those by Williams et al. [g]. It also agrees with the large 
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depinning field data of Jiang et al. (V ~ 50mV) 0. (See Fig. 8 in their paper.) We discuss 
the comparison in more detail in Sec. IX. 

The first equation holds also for the case of a sliding Wigner crystal without a magnetic 
field, in which case there is simply no Hall component. The same equations also hold when 
there is coupling to free carriers as long as we are in the large velocity regime. Figs. 1 and 
2 given previously in this paper depict the theoretical results in the Hall bar geometry. 

B). Corbino geometry 

In Corbino geometry, the electric field is controlled by virtue of symmetry of the exper- 
imental set-up. If a constant DC field E is along the x-direction, i.e., E y = 0, then for 
uj c » A, 

2 x 

Jx=—-z{E x + a dia - — -), (111) 

Pm^c Pc A 

and 

Jy = ^-(E x -2a dls ^-). (112) 

Pm U c p c U) c 

And the total current j is given to order a.di s by: 

j = ^-(E-E Th -). (113) 

Pm ^c ^c 



C). Hall Angle 

The Hall angle, i.e., the angle between the current and the external field, increases with 
the external field. In the theory in which the velocity on the right-hand-side of Eq. |31| is the 
self-consistent velocity, taking into account the disorder effects on the velocity, rather than 
the bare v in Eq. [L3|, the Hall angle becomes 45° when E/E Th = \?2u> c /(uj c — A). 

D). Shapiro anomaly 
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— * — * 

The longitudinal DC field changes by ~ E Th J n {G-VAc/^Ac) s 9 n {\G-vr,c\ — n ^Ac) across 
the resonance as the DC current or the AC frequency is swept. 

E). AC response 

The out-of-phase AC response only occurs in the longitudinal resistivity. It changes from 
inductive to capacitive as the frequency increases through resonance. 

The in-phase AC field will change by ~ E th sgn{uiAc ~ \G ■ Vdc\)\G ■ vac\/^ac) across 
the resonance. If the DC current and the AC current are not in the same direction, it will 
have both diagonal and off-diagonal components. 

B. Characteristic Lengths and Frequencies 

The elastic medium theory that we adopted to describe the sliding Wigner crystal is only 
valid in the long wavelength limit. In this subsection, we estimate various characteristic 
length scales and frequency scales that enter our previous analysis. We use numbers that 
are typical of a high mobility modulation-doped GaAs/AlGaAs heterojunction. Therefore 
m = 0.067 m e , e = 13, electron density n = 10 11 cm~ 2 , magnetic field B = 20 Tesla, 
mobility fi = 10 7 cm 2 /V ■ s, and when needed, the size of the sample is assumed to be 
3 mm x 3 mm. The ratio of the shear modulus versus the bulk modulus a is taken to 



be its classical zero-temperature value of 0.02 [pT9| . Of these parameters, a has the largest 
uncertainty, with corrections from both quantum and thermal effects, the former being 
perhaps more important near the zero-temperature Wigner solid/FQHE liquid boundary, 
the latter more important near the thermal melting transition of the Wigner solid. 

With these numbers, the cyclotron frequency u c ~ 10 13 s -1 , Q ~ 10 12 s -1 . At qa = 
1/3, the zero-field phonon frequencies are: u>l ~ 10 12 s _1 , uo^ ~ lO 11 ^^ 1 , and the acoustic 
magnetophonon frequency is ~ 10 10 s _1 . The collision frequency deduced from the zero-field 
mobility is r _1 ~ 10 9 s _1 . How this number should change in a strong magnetic field is an 
open question, as is the relation between r _1 and the parameter A in Eq. [I]. Values of A 
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inferred from the high field I — V measurements are much larger, A ~ 10 12 — 10 13 s _1 . 

Typically, measurements are done with a current / ~ InA in the sliding regime. Assum- 
ing a large fraction of the Wigner crystal is sliding, we have v ~ 10~ 2 m$ -1 . This gives a 
washboard frequency of u v between 10 5 to 10 6 s _1 . We see that this is much smaller than 
the characteristic phonon frequencies involved. Therefore, the elastic medium theory should 
be valid. We also have ^ ~ 10~ 4 , implying a typical wavelength 10 4 a, which approaches the 
typical size of a sample. It has been shown that within the perturbation theory for a CDW 
system ||12| , the broad band noise vanishes upon average over volume in an infinite system, 



but does not vanish if the system size is not infinite. The fact that k w is comparable to 
system size in the present case may be responsible for the sizable broad band noise observed 
experimentally |||| . It would be very interesting to see if experimentally the narrow band 
noise around u v could be observed. K\ is between 1CT 3 to 10~ 4 , depending on the size of A 
chosen. 

We note that if one converts the resistivity in the sliding state as measured by the 
experiments 0-0] to a scattering frequency using Eq. ^ in the high field limit, the scattering 
frequency A is rather large. In fact, it is comparable to u c , the largest frequency scale in 
the problem. It is much larger than the scattering rate deduced from the zero-field mobility, 
and much larger than that deduced from the metallic state resistivity at larger filling factors 
than v = 1/5 ||. At present, it is unclear what mechanism is responsible for such a high 
resistivity in the sliding regime. 

We suggest that it is, at least at moderately high temperatures, due to the damping 
caused by the backflow of free carriers to screen out the charge density fluctuations in a 
sliding Wigner crystal p8 |. It is clear from Eqs. |76| and [77] and the ensuing discussion that 



the free carriers lead to a longitudinal mode dissipation proportional to the free carrier 
resistivity. Although in the present problem, this does not enter the nonlinear conductivity 
to second order in the randomness, it presumably will in higher orders, as in the CDW 
problem [E3 . 
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C. Nonzero Temperatures 

The calculations reported in this paper did not include thermal effects explicitly, while at 
least in some experiments such effects are quite apparent. In this subsection, we outline 
the effects that we expect to occur at T > 0. 

We believe the four important effects occurring as T increases are: 1) an increase in 
the number of free carriers; 2) a decrease in the strength of the interaction of the Wigner 
crystal with the impurity potential; 3) an increase in the number as well as in the size of 
dislocation pairs; and 4) the subsequent melting of the Wigner crystal, presumably via a 
Kosterlitz-Thouless transition [^J. 

The first point is straightforward: presumably, creation of free carriers involves excita- 
tions over an energy gap, A g , and so the number of free carriers and thus the free carrier con- 
ductivity should have a roughly activated temperature dependence. One therefore expects 
the conductivity at E < Exh to be activated. Interestingly, in at least some experiments 
the conductivity at E > Eth is also activated, with the same gap, suggesting that the 
conductivity at E > Etk is also controlled by the free carriers. An additional point, inferred 
from the experiments in 0, is that there is one single excitation gap at temperatures below 
~ lAOmK. In other words, the low-energy charged excitations in the Wigner crystal appear 
to be well-defined, and well-separated from higher energy excitations. 

Free carriers not only provide conductivity and dissipation, they may also screen the 
impurity potential. Also thermal fluctuations of the Wigner crystal may tend to average 
out the impurity potential. For these reasons, we expect that the effective strength of the 
random potential will get weaker as the temperature is raised. 

The third effect of finite temperatures, the presence and the proliferation of dislocation 
pairs, requires mode discussion. The conventional Fukuyama-Lee-Rice |16[ theory of pinning 



and nonlinear conductivity is based on a "phase-only" model in which the amplitude flutua- 
tions of the order-parameter or the dislocations are not allowed. One important consequence 
of forbidding the formation of dislocations is the existence of a well-defined threshold field 
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below which the whole solid is pinned and above which the whole solid moves collectively. 
In a model in which dislocations are allowed, it is possible for some portions of the solid to 
move while other portions remain pinned at any given moment; and the depinning transition 
is no longer sharp nor is it unambiguously defined |14|]. In a two-dimensional ordered solid, 



the elastic energy of an isolated dislocation is logarithmically divergent in system size so at 
low temperatures the dislocations are created via thermal activation of bound pairs whereby 
the divergent elastic energy due to each one cancels out. As T increases, the number of 
pairs increases as does the typical size of a pair, C,kt- Dislocation pairs proliferate and £kt 
diverges as the Kosterlitz-Thouless melting temperature is approached from below |29| . 



Another way to see this is as follows. At zero or very low temperatures, disorder pinning 
dominates over the thermal effects and the Fukuyama-Lee-Rice theory |16|] gives the disorder- 



induced correlation length £y beyond which the positional order disappears. As temperature 
is raised, £t decreases according to £t oc c t , proportional to the (decreasing) shear modulus, 
assuming that the change in disorder potential itself is smooth and featureless. However the 
increasing thermal fluctuations also give rise to a thermally-induced correlation length £kt 
on the order of the size of a typical dislocation pair in the Wigner crystal. £xr increases very 
sharply as the melting transition is approached and diverges at the transition. Therefore 
it must become larger than the decreasing £t before the melting transition occurs, thus 



invalidating the Fukuyama-Lee-Rice description which considers only phase-fluctuations |L6 
We expect that as T increases and dislocation pairs proliferate and increase their size, the 
depinning transition will become rounded and a more appropriate physical picture of the 
system at E > would involve a combination of pinned and moving regions. The effects 
of dislocations on the nonlinear conductivity have been studied in numerical simulations by 
Shi and Berlinsky for vortex lattices in type-II superconductors at zero temperature ||15| . 

Finally, the issue of the melting of the Wigner crystal is also straightforward. The 
melting temperature may be roughly estimated from the Kosterlitz-Thouless theory and the 



theoretical shear modulus [|29|. It is found to be of order IK by Fisher [E0], much higher 



than the temperature at which the observed Etk vanishes. It is in fact much higher than 
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the temperatures at which these nonlinear transport measurements were typically carried 
out 0. Combining this with the observations that we made in the last two paragraphs, we 
conclude that it is unlikely that the vanishing of Etk observed in the experiments 0-0 was 
due to the Wigner crystal melting transition. 

IX. SUMMARY OF RECENT EXPERIMENTS AND COMPARISON TO 

THEORY 

So far we have focused on the theoretical aspects of the problem of a sliding Wigner 
crystal in a strong magnetic field when the sliding velocity is large. In what follows, we 
give a brief review of the relationship of our theoretical results to the nonlinear conductivity 
measurements reported in the literature. Seemingly similar experiments obtain apparently 
inconsistent results. We discuss these apparent inconsistencies in detail. Experiments re- 
garding the Hall resistivity were discussed earlier in this paper; so here we will focus on the 
longitudinal conductivity only. We review and compare the DC transport experiments in 
Sec. IXA, and we compare the theory and the experiments in Sec. IXB. AC + DC interfer- 
ence effects are briefly discussed in Sec. IXC. 

A. Nonlinear DC Transport Experiments 



The basic idea of these experiments is quite simple [^|-[7| JTT|JT6[| : a large enough DC electric 



field will dislodge the pinned Wigner crystal and cause it to slide. Therefore, by measuring 
the conductivity as a function of the applied DC field, one expects to observe a sudden 
increase of conductivity, and the accompanying noise generated by the disorder scattering 
off the sliding Wigner crystal. 

There are to our knowledge four sets of published results probing the nonlinear conduc- 
tivity of the insulating phases around the filling factor v = 1/5 fractional quantum Hall 
state: Goldman et al, Phys. Rev. Lett. 65, 2189 (1990) (Ref. |§); Williams et al, Phys. 
Rev. Lett. 66, 3285 (1991) (Ref. §); Li et al, Phys. Rev. Lett. 67, 1630 (1991) (Ref. §); 
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and, Jiang et al., Phys. Rev. B 44, 8107 (1991) (Ref. [0]). Among them, work by Jiang et 
al. |7|] is done at somewhat higher temperatures and covers a broader range of applied field. 
In the first half of this subsection, we shall discuss and compare the first three experiments; 
in the second half, we summarize the findings from the last one. 
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A. Sample quality 

The mobility of the samples used in Ref. [f§] was not given in the paper. However, a rising 
resistivity at v — 1/5 as the temperature is lowered is observed, although at temperatures 
reported in the paper the resistivity at v — 1/5 is a well-defined local minimum. The lowest 
temperature state at v — 1/5 is therefore insulating in these samples. What impact, if any, 
it may have on the insulating phases around v = 1/5 is not clear. 

The mobilities quoted in Ref. for various samples used range from 4 x 10 6 to 9 x 
10 6 cm 2 V~ l s~ 1 , and those quoted in Ref. H were somewhat lower ~ 1 x 10 6 Hall liquid 
state was observed. The size of the FQHE gap at v — 1/5 was reported to be ~ 0.1K in 
Ref. 0; it was not given in Ref. ||. 

B. Threshold field, broad band noise, and their temperature dependence 

The observed apparent depinning threshold field is ~ \mVcm~ 1 in Ref. [|J, which is 
similar to that in Ref. ||, found to be ~ l.2mVcm~ l , but differs greatly from the much larger 
value of ~ 250mVcm~ 1 found in Ref. ||. It is certainly worth checking if in Ref. |5j there is 
another threshold at fields similar to those in Refs. [^,0 which may have escaped observation. 
We also note here that the conventional four-terminal measurement was done in Ref. ||, 
but a two-terminal geometry was adopted in the nonlinear conductivity measurement in 
Ref. 0. In Ref. 0, both the four-terminal and the two-terminal measurements were carried 
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out in the low field (less than lOmVcrn^ 1 ) regime, and were found to give essentially the 
same results for the depinning threshold field and the nonlinear conductivity. 

In addition, at fields above the depinning threshold, broad band noise (BBN) in the 
frequency range of ~ IK Hz was observed and its power spectrum measured in Refs. [[|,^| . 
But in Ref. [Q, the external field at which the BBN was observed is greater by a factor of 
~ 70 than that at which the nonlinear conductivity set in. No explanation was given for 
this rather large discrepancy @]. In Ref. ||, the threshold field for the onset of BBN and 
that for the nonlinear conductivity were found to agree with each other. 

Despite the wide range of the observed depinning threshold field itself, its temperature 
dependence was remarkably similar in Refs. @-§f. It decreases smoothly and continuously 
as temperature rises, and becomes unobservable at about ~ 120mK. It was suggested in 
Refs. HH that the disappearance of an apparent depinning field could be due to the melting 
of the Wigner crystal at this temperature. The BBN in Refs. |4]|| had roughly the same 
temperature dependence as the depinning threshold field. In particular, it disappears at 
about the same temperature as the threshold field in Ref. ||, but in Ref. ||, it was found 
the disappear earlier, at ~ 60mK. It was therefore cautioned in Ref. |]] that it remains 
unclear what is the significance of these temperatures. 

We should also note that although Refs. f| and || reported similar threshold fields, 
the details of the depinning transition appear to be rather different: In Ref. @, a step- 
function-like increase in differential conductivity was observed at the depinning transition 
while a gradual and smooth one was found in Ref. 0] at and above the depinning field. The 
transition was the sharpest for the insulating phase at v = 0.214, slightly above v = 1/5 in 
Ref. || . In Ref. || with much larger depinning fields, the transition also appears very sharp 
for all filling factors reported, down to v ~ 0.15. 

C. Magnetic field dependence of the threshold field 

In Refs. m and 0], studies were carried out for the magnetic field, or the filling factor 
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dependence of the depinning threshold field at temperatures low enough so that the threshold 
field does not vary with temperature. This information is useful as such dependences near 
the transition from the insulating phase into the FQHE liquid phase at and immediately 
near v = 1/5 may yield information about the nature of the phase transition fPTfl . 

In Ref. J5| , the threshold field was found to decrease rapidly and monotonically approach- 
ing the FQHE liquid phase at v — 1/5 on both sides of the magnetic field. This differs from 
the observations in Ref. || in which the threshold field was found to increase rapidly on both 
sides of v = 1/5 as the (reentrant) insulating-FQHE transition was approached. For large 
enough field, or small enough filling factor, as the insulator-FQHE transition was far enough, 
the depinning threshold field was found to increase with the magnetic field in Ref. |J as was 
in Ref. §. 

D. Temperature dependence of the differential resistivity above the apparent threshold field 
In Ref. |J, only the experiments done at the lowest temperatures (~ 22mK) were re- 
ported, and it was found that the differential resistivity is T-independent above the threshold 
field at this temperature. In Ref. with a much larger depinning threshold field, a wider 
temperature range (between 20mK and lOOmK) was examined. It was found also that 
the differential resistivity is T-independent above the depinning transition, but strongly In- 
dependent below it. From these results, we may conclude that above the threshold field, 
conduction is no longer activated. We may also conclude that the main damping mechanism 
above the threshold field should be the same for all temperatures in [H. 

In Ref. [U, the differential resistivity remains field-dependent above the depinning tran- 
sition and the temperature dependence is difficult to establish clearly with the data given 
there. 

E. Size of the insulating gap 

Below the depinning threshold field, the conduction in the insulating phase is presumably 
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due to thermally excited charged defects/excitations. If there is only a single mechanism for 
the thermal excitations, the measurement of the temperature dependence of the conductivity 
yields the size of the gap for such excitations. 

The size of the insulating gap was not reported in Refs. ]5],|6[]. In Ref. ||], the tempera- 
ture dependence of the longitudinal resistivity in the insulating phases shows a complicated 
behavior, apparently involving two regimes with different temperature dependence in the 
temperature range of 30 — 100 mK. The insulating transport gap was reported to be ~ 0.5K. 
We note that this gap is much larger than the critical temperature (~ 0.12K) at which the 
threshold field disappears, postulated to be related to the melting of the Wigner crystal [|J . 
Furthermore, the insulating gap is well-defined at temperatures above the so-called melting 
temperature. This may be taken as a sign that the amplitude of the charge density does not 
vary too much in this temperature range, since the size of the insulating gap is presumably 
strongly dependent on the charge- density wave amplitude. 



A 2 Summary of findings from Ref. @ 

We believe that the samples used in |7| are probably of higher quality than those used 
in Refs. [|]-fj, because the FQHE gap at v — 1/5 in the samples of Ref. J7| was 1.1K, much 
larger than that in the others. Interestingly, the activation gap in the insulating phases is 
larger by a factor of 3 -- 4, than that in Q where the v = 1/5 state is insulating at the 
lowest temperatures, for the same filling factors. If the size of the insulating gap can be used 
as a measure of the strength of the insulating phases, we then must conclude that as sample 
quality improves, or as disorder becomes weaker, the insulating phases become stronger. 
We believe this is strong evidence supporting the notion that these insulating phases are 
primarily due to interaction, and not disorder. 

Threshold nonlinearity in conductivity was also observed in M for applied voltages 
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~ 0.1 mV across samples of typical dimensions 3mm x 3mm, giving a depinning field 
~ 1.3mV cm -1 , similar to those in || and ||. However, at and above the threshold field, 
there is no clear step- wise change in the resistivity. This is similar to |4]], but differs from 
H in which lower temperatures were reached. 

Although the threshold field diminishes at a temperature about 100 mK, below and well 
above the threshold field (within a factor of 10 of Exh)-, Pxx remains thermally activated. 
This differs qualitatively from the observations in || which are done at generally lower 
temperatures. Furthermore, the conduction activation gap remains the same below and well 
above the threshold field. In addition, not only the gap itself is an order of magnitude larger 
than the so-called critical temperature, but also the activated behavior extends above the 
critical temperature by at least a factor of 3 - - 4. Two related conclusions that one can 
draw from this important observation are: 1). The dominant damping mechanism giving 
rise to resistivity above and below the low-field threshold field ought to be the same; and 2). 
The solid retains its identity up to temperatures higher than the critical temperature by at 
least a factor of 3 - 4. The first point has a natural explanation from the coupled motion 
of a sliding Wigner crystal with thermally activated carriers. The second point supports 
the notion that the critical T at which the apparent threshold disappears does not signify 
the melting temperature of the Wigner crystal. We offer a speculation for this experimental 
observation in the next subsection. 

A voltage breakdown is observed in the high field limit around 50 mV in J7|, in the 
same range as where the work in || was carried out. It is suggested by the authors of |7j 
that the power dissipation observed in JFj is consistent with an excessive electron heating 
interpretaion. 

In addition, it is also found in |7| that the activation gap in the insulating phases may 
change upon thermal cycling to room temperature by a factor of two, as does the conductivity 
threshold. Upon thermal cycling to room temperature, the disorder configuration must 
change due to diffusion or other thermally activated processes. This suggests that disorder 
also plays an important role in determining the insulating gap. There is also a correlation 
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between the density and the activation gap, which may be taken as evidence for the role 
of interaction. The insulating gap decreases toward v = 1/5 in its vicinity, similar to the 
work in |4j] . Although the activation gap itself changes upon thermal cycling, its variation 
with filling factor is the same within a cycle. Therefore the filling-factor dependence of the 
insulating gap around the v = 1/5 fractional quantum Hall liquid state may be intrinsic to 
a clean system. 

We see from the above discussion that the experimental situation is not very clear, 
and different experiments often contradict each other in important details. In the next 
subsection, we attempt to interprete some of these experiments related to the nonlinear 
conductivity in light of our theoretical results. 

B. Comparison of Theory and Experiments in the DC limit 

According to the behavior of the conductivity at fields above threshold, the current 
experiments reviewed above fall into two categories: 

1). The differential conductivity is voltage- and temperature- independent above threshold. 
This is reported in Ref. H which finds a very large threshold field (~ 50m V) at the lowest 
temperatures ~ 36mK, and in Ref. || which finds a much smaller threshold field (~ 0.2mV) 
and is done at very low temperatures ~ 20mK. 

2). The differential conductivity above the threshold field continues to depend on voltage 
and temperature, as reported in Ref. [§J and Ref. f7j. The threshold field is found to be 
small (~ 0.2mV) [f|,[7|]. The temperature dependence of the differential resistivity is found in 
Ref. |7| to be thermally activated with an activation gap identical to that below threshold, 
for temperatures between ~ 80 mK and ~ 200 mK. However, it is also reported in 
Ref. [7] that there exists a second threshold field at ~ 30mV above which the differential 
resistivity is no longer field- or temperature-dependent, and its value is consistent with those 
in Ref. @,|]. 

We first compare our theory with experiments done at the lowest temperatures where 
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thermally activated free carriers are unimportant, i.e., the experiments in the first category 



described above. The predicted I — V curve from Eq. |109| agrees very well with the low 
temperature experiments (Compare Figs. 1 and 2 from this work to Fig. 2 in Ref. || and 
Fig. 2 in Ref. ||) - that is, the differential DC resistivity does not depend on the driving 
field, or equivalently on the resulting current above the sliding threshold. In addition, the 
intercept of the I — V curve from the high velocity perturbation theory gives also a measure 
of the sliding threshold field. From Eqs. |109| and |111| , one sees that this statement can only 
be true when, in plotting the I — V, the field V is measured along the direction of the current 
flow. 

We now consider the large field (> 30m V), therefore high velocity, data reported in 
Ref. in which electron heating was invoked to explain the data. We wish to propose 
an alternative possibility. The differential conductivity in the V > 30mV regime is field 
and temperature independent, and has roughly the same value (~ 10 5 Ohms) as the lowest 
temperature data in Ref. ||. This experimental fact |JT| is consistent with our result that 
the large- velocity differential resistivity is ^^.-independent. We propose that the second 
apparent threshold field in fact signifies a crossover from a regime where disorder effects 
enter non-perturbatively to one where they can be treated perturbatively as was done in 
this paper. As the velocity decreases, it is expected that our second-order perturbation 
theory becomes no longer valid, and that a£ x affects the differential conductivity, In this 
limit, the effective damping constant may be of order Ai oc (o - ^.) -1 , and would lead to an 
activated temperature dependence as observed. This is, however, beyond the range we can 
treat by the present perturbation theory. 

Our predicted V — I characteristic along the direction of the current j for a sliding 
Wigner crystal: E = p m Xj/pl + Etk agrees very well with the experiments in III done 
at the lowest temperatures and the high-field results in II and 0. The details of the 
pinning-depinning transition are not accessible from the present perturbation theory and 
there is no clear picture from the experiments. The strong temperature dependence of 
the differential conductivity above the small threshold field observed in |7j at somewhat 
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elevated temperatures indicates the possible role of thermally excited free carriers in the 
sliding regime beyond that revealed by our perturbation theory. 

We finally speculate on the possible cause of the disappearence of the apprarent threshold 
field at temperatures around 100 — 200 mK observed in several recent experiments @-0]. It 
was suggested that it may signify the melting of the Wigner crystal in Refs. ||||, but this 
suggestion was called into question in Refs. |||7j . We here wish to echo the cautions that have 
already been raised. First of all, one could estimate roughly the classical melting temperature 
from the shear modulus, assuming a Kosterlitz-Thouless (KT) melting mechanism |2J| . We 
find a temperature Tkt around 400 — 600 mK. Secondly, the magneto-optical measurements 
3Cfl indicate the possibility of a still higher melting temperature oi 1 — 2 K. Thirdly, above 



the so-determined "melting temperature" , the measured I — V curves continue to be highly 
nonlinear [§],[F|], which is also inconsistent with the notion that the Wigner crystal has already 
melted into a liquid. A natural alternative to this interpretation is to assume what has 
been observed is in fact the thermal depinning of the Wigner crystal. At temperatures 
above 100 — 200 mK in these experiments, we propose that thermal fluctuations destroy the 
pinning due to disorder potential, and the subsequent sliding motion is one of a solid with 
large thermal fluctuations and subjected to disorder. 

C. Theory and Experiments for AC + DC Interferences 



Finally, we briefly remark on some recent experimental work |32|j33|1 on the AC + DC 
interference effects in a sliding Wigner crystal around filling factor v — 1/5. In the absence 
of a DC current, the AC experiments by Li et al, in the frequency range of ~ 30 — 100 



MHz have established the capacitive response of the pinned Wigner crystal [^ . A detailed 
analysis and the possible implications for the nature of the insulating phases from the pure 



AC measurements were given in their original paper |32 . 

So far, there have been no reports on the observation of the Shapiro anomaly in the DC 
response due to an AC field. However, preliminary data have established the interference 
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effects in the AC response in the presence of a large DC current. This situation is not 



surprising since from the experience with the CDW systems | TI[| , the requirements on the 
sample quality are usually less stringent for the observation of the latter effect. 

In Fig. 3, we give a qualitative sketch of the calculated inductive anomaly from Eqs. [59| 
and p.01| in Sec. VIB. The dotted portion of the curve indicates that our linear AC response 
theory is not reliable at AC frequencies very close to the washboard frequency of the DC 
current. Nonetheless, its overall shape is in good agreement with the experimental observa- 
tions, where it was found that the capacitive AC response in the absence of the DC current 
becomes inductive in its presence at frequencies below 2 MHz ||33| . However, a direct tran- 
sition from the inductive response to capacitive response as the AC frequency increases has 
not been observed due to the limited AC frequency range attainable in the interference effect 
measurements |33| . 

Our theory (see Eq. ^ and Eq. |103|) also suggests that at low AC frequencies, the linear 
slope of Imp(uj) in uj at small uj would depend on the DC current as ~ 1/ j 2 DC . To observe 
this effect experimentally, care must be taken to ensure that the AC-response remains in the 
linear regime, i.e., the AC-component amplitude must be reduced along with its frequency uj 
to ensure the accuracy of the small parameter expansion of the Bessel's function (see Eq. [JB| 
and Eq. |97|). 

Clearly, further experiments are needed to test other aspects of our theory in terms of the 
AC + DC interference effects. We refer the reader to Sec. VI for a more detailed description 
of what one might expect from the present perturbation theory. 

X. CONCLUSIONS 

In conclusion, we have studied the sliding motion of a Wigner crystal in a strong magnetic 
field. We obtain the form of the nonlinear resistivity in the regime of large sliding velocity. 
The Hall resistance of a sliding Wigner crystal is found not to be changed by disorder 
scattering. We predict the AC + DC interference effects and compare the present case to 
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the CDW systems and point out the differences. We also give a brief summary of the available 
experiments that measure the nonlinear conductivity of the insulating phases around filling 
factor v — 1/5. At low temperatures and/or large depinning fields, the experimentally 
observed nonlinear I — V curves are in agreement with our theoretical results. In the 
regime of smaller fields and elevated temperatures, results from various experiments differ 
in important ways as we discussed in the paper, and a definitive comparison to the present 
theory is difficult to make. 
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Figure Captions 

Fig. 1. Theoretical nonlinear I— V properties of a sliding Wigner crystal in a strong magnetic 
field in the presence of disorder for the Hall bar geometry. / is the total current, Vl is the field 
component parallel to I, and Vh is the field component perpendicular to /. The crossing 
point of the two curves is where the Hall angle becomes 45°. The dotted portion in Vl 
indicates that the present perturbation theory is not reliable when the external field is only 
slightly greater than the depinning threshold field, indicated by Vtk- We have assumed uj c 
greater than A. The dashed portion is the conductivity due to free carriers at fields below 
the threshold field. It is extended to fields slightly above Vth f° r clarity of illustration. The 
line of Vh as a function of / is valid and extends all the way to zero / provided that the 
transport current is due to the sliding of the whole Wigner crystal under a combined total 



external field J V H + V£ greater than Vth- 

Fig. 2. The theoretical differential resistivity of a sliding Wigner crystal in a strong magnetic 
field as a function of the external field for the Hall bar geometry. The depinning threshold 
field is indicated by Vth- We have assumed u c greater than A. The dotted portion of the 
dVi/dl curve indicates that the theory is not reliable near threshold field. The dashed 
portion of the dV^jdl curve indicates the conduction of thermally excited free carriers in 
this regime where dV^jdl = p^ x , the diagonal resistivity of the free carriers. 

Fig. 3. The inductive anomaly in the imaginary part of the AC response of a Wigner crystal 
in the presence of a large DC current. The dotted portion indicates that the linear AC 
response theory is not reliable near frequency locking. 
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